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Abstract 

We investigate the superconductivity in Sr2Ru04 on the basis of the three-dimensional three- 
band Hubbard model. We propose a model with Coulomb interactions among the electrons on the 
nearest-neighbor Ru sites. In our model the intersite Coulomb repulsion and exchange coupling can 
work as the effective interaction for the spin-triplet paring. This effective interaction is enhanced 
by the band hybridization, which is mediated by the interlayer transfers. We investigate the 
possibility of this mechanism in the ground state and find that the orbital dependent spin-triplet 
superconductivity is more stable than the spin-singlet one for realistic parameters. This spin-triplet 
superconducting state has horizontal line nodes on the Fermi surface. 

PACS numbers: 74.20.Fg, 74.20.Mn, 74.20.Rp 
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I. INTRODUCTION 



The nature of the superconductivity in Sr 2 RuC>4 has drawn much attention since its 
discovery in 1994.^2 A lot of experiments have provided evidence that the superconduc- 
tivity is unconventional. For instance, the superconductivity is extremely sensitive to the 
non-magnetic impurity scattering in contrast to Anderson's theorem on a conventional su- 
perconductor.— Miyake and Narikiyo have successfully shown that such an anomalous effect 
of impurity in Sr 2 Ru04 can be explained as an evidence of the spin-triplet pairing super- 
conductivityi Nuclear magnetic resonance (NMR) measurement has revealed that the 17 
Knight shift is almost unchanged in the transition into the superconducting phase. 5 Fur- 
thermore, muon spin relaxation (/xSR) time measurement^ and polarized-neutron scattering 
study 7 clarified that in the superconducting phase the time reversal symmetry is broken. 
From these experimental evidences it is almost confirmed that the superconductivity in 
Sr 2 RuC>4 is the spin-triplet superconductivity. In the past few years, the momentum depen- 
dence of the superconducting gap function has become the central issue of this spin-triplet 
superconductor. In Sr 2 Ru04 the Fermi surface consists of three cylindrical pieces mainly 
originated from the four Ru-4<i electrons in three t 2ff orbitals . 8 i 9 i 10 Agterberg et al. insisted 
that the temperature dependences of specific heat, penetration depth, and thermal conduc- 
tivity can be explained by the orbital dependent superconductivity— Additionally, recent 
specific-heat measurement at low temperature suggests the existence of line nodes.— 

In order to determine the location of the nodes, we need the experimental results obtained 
by directional probes. In Sr 2 Ru04 the magnetothermal conductivity measurement seems the 
most powerful tool to investigate the location of the nodes . 13 > 14 Two groups have reported 
that the thermal conductivity has no notable anisotropy when the magnetic field is applied 
to the direction parallel to the conducting plane. These results are quite different from 
the result of the cuprate superconductor, and they suggest that the pairing state with 
vertical line nodes has less possibility for the candidate in Sr 2 RuC>4. Thus the paring state 
with horizontal line nodes seems to be appropriate to explain these experimental results for 
Sr 2 Ru0 4 . 

Since Sr 2 Ru04 has single-layered perovskite structure as in the case for La 2 _ x Sr x .Cu04, it 
has been supposed that its superconductivity is mediated by largely enhanced fluctuations 
common to these two-dimensional materials. 15 i 16 i 17 i 18 i 19 However, it seems difficult to explain 



2 



the spin-triplet paring state with horizontal line node. In order to solve this problem, 
Hasegawa et al. listed the possible odd-parity states on the basis of the group-theoretical 
analysis.— In their analysis they took notice of the body-centered-tetragonal lattice of Ru 
with lattice constants a and c. And they insisted that in order to stabilize the gap function 
with the horizontal line node the effective interaction for electrons at r and r ± (a/2)x ± 
(a/2)y ± (c/2)z is crucial. Zhitomirsky and Rice have successfully shown that the gap 
function with the horizontal line node may lead to the temperature dependence of the 
specific heat observed in experiments.- 21 Futhermore, Annett et al. have reproduced the 
experimental data of the superfluid density and the thermal conductivity on the basis of the 
multiband attractive Hubbard model with interlayer coupling.— 

In this paper we propose that the superconductivity in Sr 2 Ru0 4 is mediated by Coulomb 
scatterings among the electrons at r and r ± (a/2)x± (a/2)y ± (c/2)z. Our model Hamilto- 
nian is the three-dimensional (3D) three-band Hubbard model with quasi-two-dimensional 
character. Our microscopic description of the superconductivity in Sr 2 Ru04 may be con- 
sidered as an application of the two-band mechanism superconductivity to the spin-triplet 
Cooper pairing,— or as the three-dimensional version of the spin-triplet superconductivity 
in the one-dimensional chain with long-range attractive Coulomb interactions.— 

II. 3D THREE-BAND HUBBARD MODEL 

We consider three t 2g orbitals of Ru-4d electron, i.e., d xy , d yz , and d zx , in our 3D three- 
band Hubbard model. It is represented in real space as 

rr'tptp'tj 

+ X/ Upp' ( r , r ) C } pra C^^^dp'r'a' C<_p ra 
rr'ipip'acr' 

+ X] ^v'( r > r ) C ^rcr C ^'r'(T' C ¥'r'(7'C (/ p'rcr 
rr'ipip'crcr' 

+ Yl ^^( r > r ')4r<7 C trV'Vr' CT 'Vr<T ) (!) 
rr'ipip'acr' 

where c vra (c^ rcr ) is the annihilation (creation) operator of d electron with orbital tp = 
{xy,yz, zx} and spin a = {f, 1} on site r. e v are site energies, as we set £ Z x( y z) = A > 
and e xy = 0. t w '(r, r') are hopping integrals, as set 

t zxzx (r, r ± ax) = t yzyz (r, r ± ay) = t , (2) 



t Z xzx(r, r ± ay) = t yzyz (r, r ± ax) = t x , (3) 
^zx yz (yz zx) (r, r i ax + ay) 

= -t zxyz ( yzzx ){r, r ± ax - ay) = ±t 2 , (4) 

t X yxy(r, r ± ax) = t xyxy (r, r ± ay) = £ 3 , (5) 

t xyxy (r,r ±ax±ay) = t 4 , (6) 

/ a a c \ . 

****«:(»*!,*) (r,r± -x± -y ± -zj = t ± , (7) 



and 
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/ a ^ a A 

^zxxy(xyzx) ^ 2~^ 2^ 

tzxxy(xyzx) I r ; r i ~^ 2~ Z 
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tzxxy(xyzx) I T i ~X 2~^ ~^ 2~ Z 

/ a * . a „ C ^ 

= tyzxy(xyyz) I I", I" + - X ^ ^ 2~ Z 

/ a „ . a „ c „ 

= ^yzxy{xyyz) I r , r — 2~ X 2^ ~ 2~ Z 

= — ^y^xj/ (zj/i/z) ( i", r + 2 X 2^ — 2 Z 
/ a „ a „ c „ 

= ~~ ^j/zxj/ (xj/j/x) ( r ? r ~~ 2 X 2^ ~^ 2 Z 

= *'L (8) 

Hereafter, we only consider the on-site interactions and the interactions among the near- 
est neighbors along the c axis, because the interactions among the nearest neighbors 
on the conduction ab plane are negligible due to screening. If we take {fj} i=lj 8 = 
{[±(a/2)x, ±(a/2)y, ±(c/2)z]}, the Coulomb integrals in Eq. (0) turn out 
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U$(t, r') = U° W ,{1 - <W<W)<W + V x w E <W<, (9) 

i=i 

J-:(r, r') = J^,(l - <^)<W + ^ E <W„ (10) 

i=l 

KZ£,(r, r') = - V)^rr' + ^ E <W«) <W, (n) 

where U vv >, and are Coulomb repulsions, exchange interactions, and pair hop- 
pings, respectively. Then we transform our Hamiltonian from the representation in real 



space into the one in momentum k space by Fourier transform, and decompose it into 
H = Ho + H'. The noninteracting part Hq is represented by 
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In Eq. (|12|) we denote 
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(13) 
(14) 
(15) 
(16) 

(17) 

(18) 



and tii k = — 4^2 sin fc^. sin fc^, taking the in-plane lattice constant as unity. We can diagonalize 
H with respect to the band indices ( = {a, (3, 7} as H = J^ kCT X< £ c k a |* ko- a Cko- by orthogonal 
transformations, c^, kcr = J2( Rcp ka C^ anc ^ c v kCT = X< -R<Xk a Ckcr- The interacting part H' i 
represented by 



is 



H' 
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(19) 



where N is the number of k-space points in the first Brillouin zone (FBZ), and 

= W lv> cos T cos I cos 7T> ( 2 °) 
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J w >d = 87^, cos y cos -J cos— , (21) 
= ^ cos | cos | cos f. (22) 



III. SPIN-TRIPLET SUPERCONDUCTIVITY 

For our model we get a self-consistency equation for a gap function of the ( band, A^ k , 
within the weak-coupling formalism: 

A Ck = -\ e w /= Ac :r IA , 2 , (23) 

z k'c V (ec'k' ~ N + l^c'k'l 

where // is the chemical potential. Since our model does not include any asymmetrical inter- 
actions for spin state, e.g., spin-orbit interaction, this self-consistency equation is applicable 
to both spin-singlet and spin-triplet pairs in similar ways. For example, when we apply 
Eq. to a spin-triplet pair taking its odd parity, i.e., A^_ k = — A^ k , into account, we get 
the expression of V^i as below: 
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N . 

tfxpr 



+RtpC;\iR<p'{k'J<p V >k-k'Rc;'(pU.Rc'tp'\i' ■ (24) 



On the other hand, in the case for a spin-singlet pair, V^> can be expressed as 
Vcc'kk' = ttYI {^Ck-R^'Ck' iyl^ + E/^/k-k') -RcV'k-RcVk' 



N , 



+-R^k-R^ck' ^^'(1 — <W) + Ky.p'-k-k' R ( V'k-RcV'k' } • (25) 



When the gap magnitude A sc is small compared to band parameters, we can reduce Eq. 
into 

A Ck = In A sc E %'kk"5 (eck' - At) A f / k / , (26) 
k'c 

according to the Kondo's argument.-^ We choose our tight-binding band parameters as in 
Table HI where we take to as a unit of energy estimated as about leV. We choose them so that 
we can well reproduce the Fermi surface measured by the de Haas- van Alp hen effec t 2 ^ 2 ? 1 ^ 
as shown in Fig. Here we treat our tight-binding band parameters and Coulomb integrals 
as phenomenological ones. Thus it can be thought that our Fermi surface includes the band 
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TABLE I: Transfers and Coulomb interactions. 
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FIG. 1: Fermi surface in the case with A = 0.50. Band indices a, (3, and 7 are indicated here. 

renormalization effects due to the electron correlation, and that the Coulomb integrals are 
effective interactions reduced by Hartree-Fock decoupling. Hartree-Fock decoupling also 
affects on-site energies, which we can control by varying A. Our calculations are executed 
on equally spaced 256 3 k points in FBZ for each band. When we take 224 3 k points instead, 
our results of In A sc vary less than 3%. 

When we solve our reduced self-consistency equation (f2T)|) . we find that the spin-triplet 
state is more stable than spin-singlet ones. One of this reason is that [/* „/k-k' and J^/k-k' 
in Eq. (|24J) can always change their signs due to their wave-vector dependences as shown in 
eqs. (j20j) and ()21|) . Added to this, the band hybridization enhances the effective interaction of 
Eq. (}2"4"j) via the matrix elements of orthogonal transformations, -R^k and -R^k; m Eq- (|24p. 
As a result for these, U}, , k _ k , and J^/k-k' wor k like strong pair tunneling interactions 
among hybridized bands for spin-triplet pairing. As shown afterward, this hybridization is 
much important for our triplet superconductivity. One of the most stable pairing functions 
is 

A* = C c sin cos cos (27) 



7 



where Cf is real and takes different value on each ( band. Taking account of the spatial 
symmetry of our model, the other most stable function is 



Ay t — Y • ^% k z C 

\ = Cr sin — cos — cos — . 
Ck ; 2 2 2 



(28) 



It has indeed the same result of In A sc as the function, Eq. (}2Tj) . These pairing functions, 
Eqs. (J2*7|) and (|25)l. have been proposed as candidates of the most stable state by Hasegawa 
et a/.— In order to clarify the importance of the band hybridization, let us calculate the 
integrated effective matrix elements for our spin-triplet pairing function, Eq. (|27|). 



(29) 



where X)k ' denotes the momentum summation on the Fermi surface and A^ k denotes the 
normalized function of (}2Tj) determined by 



E 12 



k 

For example, when A = 0.50, we obtain 

' 1.357 x 10~ 3 

V ;,, /' ; ; r - 0.1057 

V tfy* v lf 3 t> 77 y y 0.4129 x 10~ 3 



1. 



(30) 



^ Vaa ^a"/ 



-1.175 x 10~ 3 
2.093 x 10~ 3 
0.04082 



-0.09952 N 
1.074 x 10~ 3 

-2.015 x io 3 y 



(31) 



Here, we can notice that the elements among the differnt bands v a7 , v^p, and vp a have 
larger absolute values than the others. This is caused by the pair tunneling between the 
different bands, which is enhanced by the band hybridization. If we hope to increase our 
spin-triplet pairing instability, we should use these elements effectively. Judging from the 
inequalities, v aj < < v 7 p < vp a , if C a ■ Cp < and Cp ■ C 7 < and C 7 • C a > 0, we 
expect that the eigenvalue of Eq. ([2 tip. (In A sc ) _1 , can take a large negative value for our gap 
function, Eq. (J27J). A large negative (lnA sc ) _1 results large A sc . Indeed, our numerically 
obtained solution of shown in Table |H] satisfy the above inequalities. Hence the pair 
tunneling enhanced by the band hybridization plays a significant role to realize our spin- 
triplet superconductivity. 

Hereafter, we assume that the order parameter of spin-triplet superconductor with three 
components (d vector) is parallel to the z axis, d(k) oc z(k x +ik y ) .— Then, we can reasonably 
construct our d vector as d z (k) = A^ k + iA^ k , which is a linear combination of our obtained 
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0.1625 


-0.1633 


0.08613 



TABLE II: C c of Eq. © in the case with A = 0.50. 




K x 

FIG. 2: Schematic pictures of gap amplitude on the Fermi surface of each band in the case with 
A = 0.50. The amplitude of each band is normalized in convenience. 

functions, Eqs. (|27)l and (J2BJ)- We can show that the amplitude of d vector vary as |d 2 (k)| oc 
yd — cos k x cos k y \cos(k z c/2)\, shown in Fig.EJ All of them have holizontal line nodes at k z = 
±7r/c and fourfold symmetries around the c axis, and their amplitudes are larger along [100] 
and [010] than [110]. These results are qualitatively consistent with the magnetothermal 
conductivity measurement s . 13 > 14 

Then we study the A-dependence of In A sc . This result is shown in Fig. El We show only 
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FIG. 3: A-dependence of In A ; 
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FIG. 4: Closeups of the Fermi surface projected on the plane with k z = 0. (a), (b), and (c) are 
in the cases with A = 0.44, 0.50, and 0.56, respectively. These areas are around the van Hove 
singular point as indicated in Fig. Q 

the case with 0.44 < A < 0.56 because in other cases In A sc becomes extremely small. We 
can point out that our superconducivity is reinforced only when 7 band has a large density 
of states. To make this situation clear, we magnify the part of Fermi surfaces and project 
it on the plane with k z — 0. We show this part of all Fermi surfaces with different A in 
Fig. |3J The large density of states of the 7 band can be realized when a piece of Fermi 
surface is close to the van Hove singular point (n, 0). We have earlier shown that the pair 
tunneling enhanced by the band hybridization plays a significant role for our spin-triplet 
superconductivity. Thus our spin-triplet superconductivity needs the two important factors. 
It might be rare that both of these two factors present simultaneously in real materials. We 
can expect that in Sr2Ru04 both of these two conditions are wonderfully satisfied. 

In our results A sc can get to e -4,084 ~ 16.8meV. And, when a piece of the Fermi surface 
becomes closer to the van Hove singular point (it, 0), A sc will be much larger. These results 
are too much larger than the experimental results of Sr 2 Ru04, estimated as 0.2 — 0.4meV. 
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This may be caused by too large estimations of U 1 and J 1 . However, we think that this 
is mainly caused by the weak-coupling formalism and neglected quasiparticles' lifetime. If 
the strong correlation effect decreases the lifetime, we should take into account the retar- 
dation effect and then A sc will be smaller. In Sr 2 Ru04 it is thought that the electrons 
correlate strongly with one another, and we should adopt the strong-coupling formalism for 
the quantitative estimation of A sc .— Although our quantitative estimation of A sc has these 
problems, as far as the whole electrons in Sr2Ru04 compose the Fermi liquid, our obtained 
gap symmetry cannot be replaced by the other symmetries. 

IV. CONCLUSION 

In this paper, we demonstrated that the spin-triplet pairing mediated by the intersite 
Coulomb scatterings is more stable than the spin-singlet one in our model. The gap func- 
tion has a fourfold symmetry and horizontal line nodes on the Fermi surface of each bands. 
These results appear qualitatively consistent with the experimental results. Therefore the 
interlayer Coulomb scatterings play a significant role in order to realize the spin-triplet su- 
perconductivity in Sr 2 Ru04. Judged from the results about superconducting gap magnitude, 
our superconductivity is much sensitive to the band parameters. Our superconductivity is 
unique to the electronic state in Sr 2 Ru0 4 , which has both the degenerated orbitals and the 
interlayer transfers among these different orbitals. 
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